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The evolution of chaos in active galaxy models with an oblate or a 
prolate dark halo component 
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The evolution of chaotic motion in a galactic dynamical model with a disk, a dense nucleus and a fiat biaxial dark halo 
component is investigated. Two cases are studied: (i) the case where the halo component is oblate and (ii) the case where 
a prolate halo is present. In both cases, numerical calculations show that the extent of the chaotic regions decreases 
exponentially as the scale-length of the dark halo increases. On the other hand, a linear relationship exists between the 
extent of the chaotic regions and the flatness parameter of the halo component. A linear relationship between the critical 
value of the angular momentum and the flatness parameter is also found. Some theoretical arguments to support the 
numerical outcomes are presented. An estimation of the degree of chaos is made by computing the Lyapunov Characteristic 
Exponents. Comparison with earlier work is also made. 

© 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinhcim 



1 Introduction 

Observational data show that disk galaxies are often sur- 
rounded by massive and extended dark matter haloes. The 
best tool to study dark matter haloes in galaxies are the rota- 
tion curves derived from neutral hydrogen (see e.g. Meurer 
et al. 1996; Dinshaw et al. 1998; Shull et al. 1998; Oppen- 
heimer et al. 2001; McLin et al. 2002; Penton et al. 2002; 
Steidel et al. 2002; Cote et al. 2005). On the other hand, nu- 
merical simulations suggest that dark matter galactic haloes 
are not only spherical but also may be oblate, prolate or tri- 
axial (Menitt & Fridman 1996; Cooray 2000; Kunihito et al. 
2000; Oiling & Merrifield 2000; Jing & Suto 2002; Wech- 
sler et al. 2002; Kasun & Evrard 2005; Allgood et al. 2006; 
Capuzzo-Dolcetta et al. 2007; Wang et al. 2009; Evans et 
al. 2009; Caranicolas & Zotos 2009). The variety of the 
shapes of galactic haloes strongly indicates that the struc- 
ture of these objects plays an important role in the orbital 
behavior and, generally, in the dynamics of a galaxy. 

In two earlier papers (Caranicolas 1997; Papadopoulos 
& Caranicolas 2006) we have studied axially symmetric or 
non axially symmetric active galaxy models with an addi- 
tional spherical halo component. In both cases, it was ob- 
served that the presence of the spherical halo had as a re- 
sult to reduce the area in the phase space occupied by the 
chaotic orbits. Therefore, it would be of significant interest 
to investigate the behavior of orbits in an active galaxy with 
a biaxial halo. On this basis, we have decided to study the 
motion in a disk galaxy model with an oblate or a prolate 
dark halo component. Particular interest will be given to the 
study of the regular or chaotic character of orbits and its 



connection to the physical parameters of the system, such 
as the flatness, the core radius of the halo component and 
the conserved component of the angular momentum. Fur- 
thermore, we shall use the Lyapunov Characteristic Expo- 
nent (L.C.E) (see Lichtenberg & Lieberman 1992) in order 
to estimate and compare the degree of chaos in each case. 

The dynamical model is presented in Sect. 2. Two cases 
are distinguished. The case when an oblate and the case 
when a prolate halo component is present. In Sect. 3 we 
study the behavior of orbits when an oblate or a prolate 
halo component is present. In the same section, numeri- 
cally found relationships between the extent of the chaotic 
regions and the physical parameters of the system are pre- 
sented. In Sect. 4 some semi-theoretical arguments are used 
in order to explain the numerically obtained results. We 
close with Sect. 5 where a discussion and the conclusions 
of this research are given and a comparison with previous 
work is made. 



2 Description of the dynamical model 

There are three components in our dynamical model. The 
disk-halo, the dense nucleus and the dark halo component. 
The disk-halo component is represented by the potential 



V d (r,z) 



-Md 
R 



(1) 
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where R 2 = b 2 + r 2 + (a + \Jh 2 + z 2 ) . Here (r, z) are 
the usual cylindrical coordinates, Md is the mass, b is the 
core radius of the disk-halo, a is the disks scale length and 
h corresponds to the disk's scale height. The dense nucleus 
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is represented by the spherical potential 

V n (r,z)= /2 , (2) 

[r z + z £ + cfj 

where M n is the mass and c„ is the scale length of the nu- 
cleus. For the dark halo component we use the logarithmic 
potential 

V h (r,z) = f\n(r 2 +(3z 2 + cl), (3) 

where f3 is the flatness parameter, while Ch stands for the 
scale length of the dark halo component. The parameter vo 
is used for the consistency of the galactic units. We have 
chosen this potential because Eqs. (1) and (2) represent a 
modified version of the analytical potential used success- 
fully by Caranicolas & Innanen (1991) to describe the mo- 
tion in an active disk galaxy. Furthermore, Eq. (3) is used 
because we believe that potential (3) is suitable for the de- 
scription of the motion in a dark halo as it produces a flat 
rotation curve (see Binney & Tremaine 2008). In this re- 
search, we shall use a system of galactic units, where the 
unit of length is 1 kpc, the unit of time is 0.977 x 10 8 yr 
and the unit of mass is 2.325 x 10 7 M Q . The velocity unit 
is 10 km/s, while G is equal to unity. In the above units we 
use the values: v = 20, a = 3, b = 6, h = 0.2, Md = 
5000, M n = 400, c„ = 0.25, while /3 and c h are treated as 
parameters. The total potential responsible for the motion of 
a test particle (star) of unit mass in the galaxy is 

V t (r,z) = V d + V n + V h . (4) 

As the total potential V t — Vt(r, z) is axially symmetric and 
the L z component of the angular momentum is conserved 
we use the effective potential 



Veff(r,z) = ^ + V t (r,z), 



(5) 



in order to study the motion in the meridian (r, z) plane. 
The equations of motion are 

9V eff ■ _ dV ef f 



r=p r , 



Pz, Pr 



dr 



Pz 



dz 



(6) 

where the dot indicates derivative with respect to the time. 
The corresponding Hamiltonian is written as 
1 



H = -i P 2 r +p 2 z ) + V eff (r,z)=E, 



(7) 



where p r and p z are the momenta, per unit mass, conjugate 
to r and z, while E is the numerical value of the Hamilto- 
nian, which is conserved. Eq. (7) is an integral of motion, 
which indicates that the total energy of the test particle is 
conserved. 

Orbit calculations are based on the numerical integra- 
tion of the equations of motion, which was made using a 
sharp Bulirsh-Stoer FORTRAN routine in double precision. 
By the term sharp we mean that all subroutines of the nu- 
merical integration code are in double precision. The accu- 
racy of the calculations was checked by the constancy of 
the energy integral, which was conserved up to the twelfth 
significant figure. 



3 Orbit calculations when a biaxial halo 
component is present 

In this section we shall study the behavior of orbits when 
a biaxial halo is present. We shall use the classical method 
of the Poincare (r,p r ), z = 0,p z > phase plane in order 
to determine the regular or chaotic character of motion. Two 
cases will be studied: (a) the system has a spherical or oblate 
halo, that is when 1 < (3 < 2 and (b) the system has a 
prolate halo, that is when 0.1 < f3 < 1. The value of Ch is 
in the range 8.5 < Ch < 21. 

3.1 Model with an oblate halo 

Figure la-b shows the (r,p r ) phase plane when (a) f3 = 1.3 
and (b) j3 = 1.8. The value of c h is 8.5, while the values of 
all the other parameters are l>o = 20, a — 3,6 = 6,h = 
0.2, M d = 5000, M n = 400, c„ = 0.25 and L z = 10. 
Here, we must note that all the initial conditions are taken 
inside the limiting curve. This curve contains all the invari- 
ant curves in the (r, p r ), z — 0, p z > phase plane and can 
be obtained from (7) by setting z — p z = (see Papadopou- 
los & Caranicolas 2006). Thus, we take the values of r and 
p r inside the limiting curve, while the value of p z is found 
from the energy integral (7). As one can see the majority 
of the phase plane is covered by chaotic orbits. Regular or- 
bits are confined mainly near the central parts of the phase 
plane. There are also some smaller islands produced by sec- 
ondary resonances. One observes that in both cases there is 
a single chaotic sea. This chaotic sea is larger in the case 
when (3 = 1.8. This suggests that the flatness parameter of 
the halo plays an important role on the character of motion. 
Our numerical calculations indicate that the chaotic region 
in disk galaxies with dense nuclei increases when the flat- 
ness parameter (3 of the halo increases, provided that all the 
other parameters are kept constant. We shall come to this 
point later in the next section. 

Figure 2a-b shows the (r,p r ) phase plane when (a) 
Ch = 11 and (b) Ch = 18.5. The value of f3 is 1.7, while 
all the other parameters are as in Fig. 1. Again, the ma- 
jority of orbits are chaotic in both cases. Here, we see that 
when Ch — 11 the chaotic sea is larger. This indicates that 
the scale length of the halo plays also an important role on 
the character of motion. We observe that the extent of the 
chaotic regions in the phase plane decreases as Ch increases, 
when all the other parameters are kept constant. Therefore, 
we conclude that the chaotic regions in active disk galax- 
ies are larger when a dense halo component is present. Fig- 
ure 3 shows the relationship between the percentage of the 
area A% covered by chaotic orbits in the phase plane and Ch 
when f3 = 1.7. We see that A% decreases exponentially as 
Ch increases. Here, we must notice that the whole area of the 
phase plane is reduced and becomes smaller as Ch increases. 

Figure 4a-d shows four representative orbits. Figure 4a 
shows a regular orbit when [3 = 1 and Ch = 8.5. The initial 
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Fig. 1 (a-b): The (r,p r ) phase plane when (a, left) ft = 1.3 and (b, right) ft = 1.8. The value of Ch is 8.5, while the 
values of all the other parameters are: v = 20, a = 3, b = 6, h = 0.2, M d = 5000, M n = 400, c n = 0.25 and L z = 10. 



40 




Fig.2 (a-b): The (r,p r ) phase plane when (a, left) Ch = 11 and (b, right) Ch = 18.5. The value of ft is 1.7, while the 
values of all the other parameters are as in Fig. 1. 




8 10 12 14 16 18 20 22 



Fig. 3 A plot of the percentage of the area A% in the 
phase plane covered by chaotic orbits vs. c/, when ft = 1.7 



The initial conditions are: r = 5,z — 0,p r o — 20. All 
orbits were calculated for a time period of 100 time units. 
The value of energy is E = 600, while the values of all 
the other parameters are: v — 20, a — 3, ft = 6, h — 
0.2, M d = 5000, M n = 400 and c„ = 0.25. The value of 
L z is equal to 10. 

In order to have an estimation of the degree of chaos 
from a different point of view, we computed the L.C.E in 
the chaotic sea for each of the corresponding (r,p r ) phase 
planes for a time period of 10 5 time units for different values 
of ft and Ch- Our numerical experiments suggest that the 
L.C.E is in the range from about 0.36 to about 0.43, when 
Ch = 8.5 and 1 < ft < 2, while it is in the range from 
about 0.40 to about 0.55, when ft = 1.7 and 8.5 < c h < 21. 
Figure 5 shows the L.C.E for the chaotic orbit of Fig. 4d. 



conditions are: r = 3, z = 0,p r0 = 30. The value of p z0 
is found from the energy integral (7) for all orbits. Figure 4b 
shows a regular orbit orbit when ft = 1.7 and c/, = 13. The 
initial conditions are: r = 6, z = 0,p r o = 0. In Fig. 4c a 
quasi periodic orbit is shown. Here, ft = 1.9 and Ch = 8.5. 
The initial conditions are: r = 9,z = 0,p r o = 0. This 
orbit is characteristic of the 4:3 resonance. A chaotic orbit 
is given in Fig. 4d. The value of ft is 1.6, while Ch — 8.5. 



3.2 Model with a prolate halo 

Figure 6a-b shows the (r,p r ) phase plane when (a) ft = 0.6 
and (b) ft = 0.2. The value of Ch is 8.5, while the values of 
all the other parameters are: vq = 20, a = 3, ft = 6, h = 
0.2, M d = 5000, M n = 400, c n = 0.25 and L z = 10. Let 
us start from Fig. 6a. In this case we see that the majority 
of the phase plane is covered by regular orbits. A consider- 
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Fig. 4 (a-d): Orbits when the system has an oblate dark halo, (a, upper left): a regular orbit when /? = 1.0 and c h = 8.5. 
The initial conditions are: r — 3, z = 0,p r o = 30. (b, upper right): a regular orbit orbit when (3 = 1.7 and Ch = 13. 
The initial conditions are: ro = 6, zq = 0,p r o = 0. (c, lower left): a quasi periodic orbit when j3 = 1.9 and Ch = 8.5. The 
initial conditions are: r = 9, z = 0,p r -o = 0. (d, /ower right): a chaotic orbit. The value of (3 is 1.6, while Ch = 8.5. The 
initial conditions are: r = 5, z = 0,p r o = 20. The value of p z0 is found from the energy integral (7) for all orbits. The 
values of all the other parameters are given in text. 




2 4 6 8 10 024G8 10 



Fig. 6 (a-b): The (r,p r ) phase plane when (a, left) (3 = 0.6 and (b, right) (3 = 0.2. The value of Ch is 8.5, while the 
values of all the other parameters are: v = 20, a = 3, (3 = 6, h = 0.2, M d = 5000, M n = 400, c„ = 0.25 and L z = 10. 



able chaotic layer is present and it is confined in the outer 
parts of the phase plane. In Fig. 6b we have the case where 
j3 = 0.2. Here again the majority of the orbits are regular 
and secondary resonances are also present. The main differ- 
ence from Fig. 6a is that now the chaotic layer has become a 
chaotic sea. This suggests that the flatness parameter of the 
prolate halo plays also an important role on the character of 
motion. We shall come to this point later in the next section. 



Figure 7a-b shows the (r,p r ) phase plane when (a) 
Ch = 11 and (b) Ch — 18.5. The value of (3 is 0.6, while 
all the other parameters are as in Fig. 1. Again, the majority 
of orbits are regular. A careful observation shows that when 
Ch = 11 the chaotic sea is larger than when Ch = 18.5. Note 
again that the whole area of the phase plane is reduced and 
becomes smaller as c/j increases. Figure 8 shows the rela- 
tionship between the percentage of the area A% covered by 
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Fig.7 (a-b): The (r,p r ) phase plane when (a, left) Ch = 11 and (b, right) ct = 18.5. The value of (3 is 0.6, while the 
values of all the other parameters are as in Fig. 6. 
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Fig. 5 Evolution of the L.C.E with the time for the chaotic Fig. 10 Evolution of the L.C.E with the time for the 
orbit shown in Fig. 4d. chaotic orbit shown in Fig. 9d. 
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Fig. 8 A plot of the percentage of the area A% in the 
phase plane covered by chaotic orbits vs. ct when (3 = 0.7. 
The values of all the other parameters are given in text. 

chaotic orbits in the phase plane and Ch, when f3 = 0.7. We 
observe that A% decreases exponentially as Ch increases. 
This behavior is similar to that found in Fig. 3 for the oblate 
halo component. 

In Figure 9a-d we present four representative orbits. Fig- 
ure 9a shows an orbit when (3 = 0.2 and c h = 8.5. The 
initial conditions are: r = 1.3, z = 0,p r0 = 6. This or- 



bit produces the upper small island shown in Fig. 5b. Figure 
6b shows a regular orbit when f3 = 0.3 and Ch = 8.5. The 
initial conditions are: ro = 5, zq = 0,p r o = 0. The or- 
bit is characteristic of the 2:1 resonance. In Fig. 9c a quasi 
periodic orbit is shown. Here, (3 — 0.6 and c h = 16. The 
initial conditions are: ro = 6.5, zq = 0,p r o = 0. The or- 
bit is characteristic of the 4:3 resonance. A chaotic orbit is 
given in Fig. 9d. The value of (3 is 0.2, while c h = 8.5. 
The initial conditions are: ro = 2, zq = 0,p r o = 0. All or- 
bits were calculated for a time period of 100 time units. The 
value of energy is E = 600, while the values of all the other 
parameters are: vq — 20, a = 3, /3 = 6, h = 0.2, = 
5000, M„ = 400, c„ = 0.25 and L z = 10. 

Computation of the L.C.E in the case of the prolate halo 
for different values of (3 and Ch, shows that the L.C.E is in 
the range from about 0.40 to about 0.55, when c/j = 8.5 and 
0-1 < (3 < 0.9, while it is in the range from about 0.51 to 
about 0.57, when /3 = 0.6 and 8.5 < c h < 21. Figure 10 
shows the L.C.E for the chaotic orbit of Fig. 9d. 

Let us now come to see how the percentage of the 
chaotic regions in the phase plane is connected with the 
flatness parameter (3. The results are given in Fig. 1 1 when 
0.1 < (3 < 2, that is, for both the prolate and the oblate 
halo together. There are tree different linear parts in this fig- 
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Fig. 9 (a-d): Orbits when the system has a prolate halo, (a, upper left): an orbit when /3 = 0.2 and c h = 8.5. The initial 
conditions are: r = 1.3, Zo = 0,p r o = 6. (b, upper right): a regular orbit when (3 = 0.3 and Ch — 8.5. The initial 
conditions are: ro = 5, zq = 0,p r o = 0. (c, lower left): a quasi periodic orbit is shown. Here f3 — 0.6 and Ch = 16. The 
initial conditions are: r = 6.5, z = 0,p r o = 0. (d, lower right): a chaotic orbit when f3 = 0.2 and Ch = 8.5. The initial 
conditions are: r = 2, z = 0,p r o = 0. The values of all the other parameters are given in text. 




Fig. 11 A plot of the percentage of the area A% in the 
phase plane covered by chaotic orbits vs. /?. The value of 
Ch is 8.5, while the values of all the other parameters are: 

v Q = 20, a = 3,6 = 6,h = 0.2, M d = 5000, M n = 
400, c = 0.25 and L z = 10. 



Before closing this section, we would like to present 
a relationship connecting the critical value of the angular 
momentum L zc (that is the maximum value of the angular 
momentum, for which, stars moving near the galactic plane 
are scattered to the halo displaying chaotic motion, for a 
given value of fi) and the flattening parameter /3, when all 
the other parameters are kept constant. The results, which 
were found numerically, are given in Fig. 12. The value of 
energy is E — 600, while the values of all the other pa- 
rameters are: vq = 20, a = 3,/3 = 6, h = 0.2, Md = 
5000, M n = 400, c„ = 0.25 and c h = 8.5. We see that the 
relationship between j3 and L zc is linear. Orbits starting in 
the upper part of the (/3, L zc ) plane are regular, while orbits 
starting in the lower part of this plane including the line are 
chaotic. An explanation of this behavior will be given in the 
next section. 



4 Some semi-theoretical arguments 



ure. In the first part, where 0.1 < /? < 0.9.4% decreases, 
in the second part where 0.9 < (3 < 1.0A% increases very 
rapidly, while in the third part where 1 < < 2A% also in- 
creases. An explanation for this behavior of the system will 
be given in the next section. 



In this section we shall present some theoretical arguments 
together with elementary numerical calculations in order to 
explain the numerically found relationships given in Figs. 
11 and Fig. 12. The contribution of the halo component to 
the chaotic regions observed in the (r,p r ) phase plane of 
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Fig. 12 A plot of the L zc vs. (3 when c h = 8.5. The value Fig. 13 A plot of \F zh | and e vs. {3. The values of all the 
of Ch is 8.5, while the values of all the other parameters are: other parameters are given in the text. 



v a = 20, a = 3, b = 6, h 
and c„ = 0.25. 



0.2, M d = 5000, M n 



400 



the system comes from two parts. The first part is the F z h 
force, which is the vertical force of the halo and the second 
part comes from the asymmetry of the dark halo component. 
The F z h force is 



F zh = 



+ (3z 2 + c{ 



(8) 



while the asymmetry of the dark halo component can be ex- 
pressed using the ellipticity (see Binney & Tremaine 2008). 
The ellipticity of the halo is defined as 



1 



b 

a 



(9) 



where a and b are the major and the minor axis of the biaxial 
halo respectively. For the prolate halo we find e = 1 — y/Jii, 
while for the oblate halo e = 1 — Figure 13 shows 

a plot of the \F z h\ vs. j3. As the scattering occurs near the 
nucleus it must be r < 1 and z < 1. The particular values 
of r and z are irrelevant. Here we choose r = ro = 0.2 and 
z = zo = 0.1. We see that \F z h\ increases linearly with (3. 
In the same figure a plot of e vs. (3 is given. The value of Vq 
is 20, while c h = 8.5. We observe that when 0.1 < j3 < 0.9 
we have an exponential decrease of e, while at the same 
time we have a linear increase of the \F zh \. The result of 
these two contrary actions is the observed linear decrease of 
chaos shown in Fig. 11. When 0.9 < j3 < 1.0 the decrease 
of the ellipticity can be considered linear. But at the same 
time, we have a considerable linear increase of the \F z h\, 
which has obtained larger values. The outcome is a rapid 
increase of the chaotic region. When 1 < j3 < 2, e increases 
slowly, while the dominant role is played by the \F z h\. As a 
consequence we see that the chaotic region increase to high 
values up to about A% = 70, when (3 reaches 2. 

The exponential decrease of the percentage of the area 
A% covered by chaotic orbits in the phase plane can be eas- 
ily explained by plotting the \F z h\ force, given by Eq.(8)vs. 
Ch, while keeping all the other parameters fixed. The results 
are shown in Fig. 14, where we have taken v = 20, (3 = 



Fig. 14 A plot of \F z h\ vs. Ch- The values of all the other 
parameters are given in the text. 

1.8, r = r = 0.2 and z = z = 0.1. We see that F zh 
decreases exponentially as Ch increases. 

In order to explain the results shown in Fig. 12 we use 
an analysis similar to that used in Caranicolas & Innanen 
(1991). When the star approaches the dense nucleus its mo- 
mentum in the z direction changes according to the equation 



mAv z = (F zt )At, 



(10) 



where m is the mass of the star, At is the duration of the 
encounter and (F zt ) is the total average F z force. It was 
observed that the stars deflection into higher z proceeds in 
each case cumulatively, a little more with each successive 
pass by the nucleus and not with a singe "tragic" encounter. 
It is assumed, that the star is scattered off the galactic plane 
after n > 1 encounters when the total change in the momen- 
tum in the z direction is of order of the tangential velocity 
W0 = L zc /r. Thus, we have 



m 



^Aw„ w (F zt )J2&U- 



(11) 



i=0 



i=0 



If we set to = 1, Yh=o Au « = L zc/r, Yh=o A ^ = T c and 
Pv 2 \z\ 



F zt w h 



r 2 + f3z 2 + c\ ' 



(12) 
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where k\ stands for the constant at a given point and for 
fixed values of the involved parameters and the vertical 
force coming from the nucleus and the disk components, 
in Eq. (11) we find 



fci + 



l3vj\z\T c 
+ pz 2 + c\ 



(13) 



The star must go close to the nucleus in order to be scattered. 
In this case we may set r = r = \z\ = c < 1 in (13) and 
obtain 



k + \fl, 



(14) 



where k and A are constants. Relation (14) explains the nu- 
merically found relationship between j3 and L zc . The dif- 
ferent slope in Fig. 12 for the prolate and the oblate halo 
can be explained because the values of c and T c are slightly 
different for the two types of the dark halo component. 

5 Discussion 

In this article, we have tried to study the regular or chaotic 
character of motion in an active galaxy model with a biaxial 
dark halo component. It is well known from our previous 
work in disk galaxies with dense nuclei (see Caranicolas 

6 Innanen 1991; Caranicolas & Papadopoulos 2003) that 
low angular momentum stars moving near the nucleus are 
scattered off the galactic plane displaying chaotic motion. 
Furthermore, this procedure is one of the main mechanism 
that produces chaos in galaxies (see Grosb0l 2002). 

On the other hand, earlier work (Caranicolas 1997; Pa- 
padopoulos & Caranicolas 2006) indicates that the role of a 
spherically symmetric dark halo in galaxies with dense mas- 
sive nuclei is to reduce the extent of the chaotic regions. It 
was this motive that drive us to consider a model of an active 
galaxy with a biaxial halo and study the dynamical effects 
of the additional component in the behavior of the orbits. In 
order to keep things simple, we have kept all the parame- 
ters of the model constant and studied the behavior of orbits 
varying only two basic parameters, that is the flattening pa- 
rameter j3 and the scale length Ch of the halo. In some cases 
we had to find the critical value of the angular momentum 
in order to connect it with the flattening parameter (3. 

It was found that when a biaxial halo component is 
present there is a linear relationship between A% and the 
flattening parameter f). On the contrary, the relation between 
A% and Ch is not linear but exponential. In both cases the 
numerically found results were explained using some semi- 
theoretical arguments. In the same sense we have explained 
the numerically obtained relationship between L zc and f3. 

Computation of the L.C.E shows that the degree of 
chaos is similar to that found for 3D time-dependent axi- 
ally symmetric galactic potentials (see Caranicolas & Pa- 
padopoulos 2003), where the L.C.E was about 0.5, while 
in 2D non axially symmetric potentials (Papadopoulos & 
Caranicolas 2006) the L.C.E was found larger, about 1. 



The main conclusion of this research is that the presence 
of a flattened dark halo component in an active disk galaxy 
has as a result to increase the extent of the chaotic regions 
observed in the (r,p r ) phase plane. This result is different 
from the results obtained in earlier papers, where the pres- 
ence of the spherical halo had as a result the decrease of 
the chaotic region. It is in our plans, to study the behavior 
of orbits in active galaxy models, when a triaxial dark halo 
component will be present, in the near future. 
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